SPITZER IRS spectra of Virgo early type galaxies: 
detection of stellar silicate emission. 
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ABSTRACT 

We present high signal to noise ratio Spitzer Infrared Spectrograph observations of 17 Virgo 
early- type galaxies. The galaxies were selected from those that define the colour-magnitude 
relation of the cluster, with the aim of detecting the silicate emission of their dusty, mass- 
losing evolved stars. To flux calibrate these extended sources we have devised a new procedure 
that allows us to obtain the intrinsic spectral energy distribution and to disentangle resolved 
and unresolved emission within the same object. We have found that thirteen objects of the 
sample (76%) are passively evolving galaxies with a pronounced broad silicate feature which is 
spatially extended and likely of stellar origin, in agreement with model predictions. The other 
4 objects (24%) are characterized by different levels of activity. In NGC 4486 (M 87) the line 
emission and the broad silicate emission are evidently unresolved and, given also the typical 
shape of the continuum, they likely originate in the nuclear torus. NGC 4636 shows emission 
lines superimposed on extended (i.e. stellar) silicate emission, thus pushing the percentage of 
galaxies with silicate emission to 82%. Finally, NGC 4550 and NGC 4435 are characterized 
by polycyclic aromatic hydrocarbon (PAH) and line emission, arising from a central unresolved 
region. A more detailed analysis of our sample, with updated models, will be presented in a 
forthcoming paper. 
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Bressan, Granato & Silva (1998) have sug- 
gested that the presence of dusty circumstellar 
envelopes around asymptotic giant branch(AGB) 
stars should leave a signature, a clear excess at 
10 /im, in the mid infrared (MIR) spectral region 
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of passively evolving stellar systems. Early de- 
tections of such an excess were suspected in M32 
(Impey et al. 1986) from ground based observa- 
tions, and in a few ellipticals observed with ISO- 
CAM (Bregman et al. 1998). The first unam- 
biguous confirmation of the existence of this fea- 
ture, though barely resolved, was found in the ISO 
CVF spectrum of NGC 1399 (Bressan et al. 2001). 
Since AGB stars are luminous tracers of interme- 
diate age and old stellar populations, an accurate 
analysis of this feature has been suggested as a 
complementary way to disentangle age and metal- 
licity effects among early type galaxies (Bressan et 
al. 1998; 2001). More specifically, Bressan et al.'s 
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models show that a degeneracy between metalhc- 
ity and age persists even in the MIR, since both 
age and metalUcity affect mass-loss and evolution- 
ary lifetimes on the AGB. While in the optical age 
and metallicity need to be anti-correlated to main- 
tain a feature unchanged (either colour or narrow 
band index), in the MIR it is the opposite: the 
larger dust-mass loss of a higher metallicity sim- 
ple stellar population (SSP) must be balanced by 
its older age. Thus a detailed comparison of the 
MIR and optical spectra of passively evolving sys- 
tems constitutes perhaps one of the cleanest ways 
to remove the degeneracy. Besides this simple mo- 
tivation and all other aspects connected with the 
detection of evolved mass-losing stars in passive 
systems (e.g. Athey et al. 2002), a deep look 
into the mid infrared region may reveal even tiny 
amounts of activity. In this letter we present the 
detection of extended silicate features in a sample 
of Virgo cluster early type galaxies, observed with 
the IRS instrument^ (Houck et al. 2004) of the 
Spitzer Space Telescope. (Werner et al. 2004). 

2. Observations and data reduction 

Standard Staring mode short SLl (7.5^m- 
15.3/im), SL2 (5/im-7.6/im) and long LL2 (14.1/im- 
21.3Mm), low resolution (R~64-128) IRS spectral 
observations of 17 early type galaxies, were ob- 
tained during the first Spitzer General Observer 



The IRS was a collaborative venture between Cornell Uni- 
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Table 1: Virgo galaxies observed with IRS 



Name 


Vt 


Date 


SLl/2 
60s 


LL2 
120s 


S/N 
6/^m. 


NGC 4339 


11.40 


.Jun 


06 


2005 


20 


14 


39 


NGC 4365 


9.62 


.Jan 


10 


2005 


3 


3 


57 


NGC 4371 


10.79 


Jun 


01 


2005 


9 


10 


40 


NGC 4377 


11.88 


Jun 


01 


2005 


12 


8 


54 


NGC 4382 


9.09 


Jul 


07 


2005 


3 


3 


59 


NGC 4435 


10.66 


Jun 


01 


2005 


3 


5 


35 


NGC 4442 


10.30 


Jan 


10 


2005 


3 


3 


46 


NGC 4473 


10.06 


Jun 


01 


2005 


3 


3 


55 


NGC 4474 


11.50 


Jun 


01 


2005 


20 


14 


38 


NGC 4486 


8.62 


Jun 


03 


2005 


3 


3 


80 


NGC 4550 


11.50 


Jun 


03 


2005 


20 


14 


42 


NGC 4551 


11.86 


Jun 


03 


2005 


20 


14 


47 


NGC 4564 


11.12 


Jun 


07 


2005 


4 


6 


51 


NGC 4570 


10.90 


Jun 


06 


2005 


3 


5 


42 


NGC 4621 


9.81 


Jan 


12 


2005 


3 


3 


63 


NGC 4636 


9.49 


Jul 


08 


2005 


3 


5 


30 


NGC 4660 


11.11 


Jan 


11 


2005 


3 


5 


40 



Cycle. The galaxies were selected among those 
that define the colour magnitude relation of Virgo 
cluster (Bower, Lucy & Ellis 1992). The observing 
log is given in Table 1. We also report, in columns 
4 and 5, the number of cycles of 60 and 120 sec- 
onds exposures performed with SLl/2 and LL2, 
respectively. The spectra were extracted within 
a fixed aperture (3'.'6xl8" for SL) and calibrated 
using custom made software, tested against the 
SMART software package (Higdon et al. 2004). 

The on-target exposures in each SL segment 
(e.g. SLl) also provide ~80" offset sky spectra 
in the complementary module (e.g. SL2) that 
were used to remove the sky background from 
the source spectrum in the corresponding seg- 
ment. Since for the LL module wc have ob- 
tained only LL2 observations, LL2 spectra were 
sky-subtracted by differencing observations in the 
two nod positions. 

2.1. Flux calibration for extended sources 

The IRS pipeline version S12 (and older ver- 
sions) is designed for point source flux extraction. 
We present here an alternative procedure that ex- 
ploits the large degree of symmetry that charac- 
terizes the light distribution in early type galaxies. 

We first obtained the real e~/sec to Jy conver- 
sion following the procedure outlined by Kennicutt 
et al. (2003) . We have corrected the conversion ta- 
ble provided for point sources by applying the cor- 
rections for aperture losses (ALCF) and slit losses 
(SLCF). The ALCF is due to the residual flux 
falling outside the aperture selected in the stan- 
dard calibration pipeline. To estimate the ALCF 
we used 4 calibration stars (HR 2194, HR 6606, 
HR 7341 and HR 7891) observed with Spitzer IRS. 
Using the Spitzer Science Center SPICE software 
package we have evaluated the correction result- 
ing from the average ratio of the fluxes extracted 
within the standard aperture and within twice the 
standard aperture. The SLCF correction is ap- 
plied to retrieve the real flux of an observed point 
source that hits the slit, accounting for the slit 
losses due to the point spread function of the opti- 
cal combination (SLCF). It is defined as the wave- 
length dependent ratio between the whole flux of 
a point source on the field of view and the flux 
selected by the slit to hit the detector. To ob- 
tain this correction we have simulated the point 
spread function of the system (PSF) using the 
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Spitzer-adaptcd Tiny Tim code and adopting a 
"hat" beam transmission function of the sht. 

After the ALCF and SLCF corrections were 
applied we obtained the flux received by the slit 
within a given aperture. The estimate of the flux 
emitted by an extended source within the selected 
angular aperture of the slit involves the deconvo- 
lution of the received flux with the PSF of the 
instrument. This correction is important to ob- 
tain the shape of the intrinsic spectral energy dis- 
tribution (SED) of the galaxy, because from the 
SLCF we have estimated that for a point source 
the losses due to the PSF amount to about 20% 
at 5/im and to about 40% at 15/um. Conversely, 
a uniform source will not suffer net losses. In or- 
der to recover the intrinsic SED we have convolved 
a surface brightness profile model with the PSF, 
and we have simulated the corresponding observed 
linear profile along the slits, taking into account 
the relative position angles of the slits and the 
galaxy. The adopted profile is a wavelength depen- 
dent two dimensional modified King's law (Elson 
at al. 1987): 



I = Io/ 



1 



y2 



1-7/2 



(Rc X fe/a)2 



(1) 



where X and Y are the coordinates along the ma- 
jor and minor axis of the galaxies, b/a is the ax- 
ial ratio taken from the literature. Iq, Rq and 7 
are free parameters that are functions of the wave- 
length and are obtained by fitting the observations 
with the simulated profile. In order to get an ac- 
curate determination of the parameters of the pro- 
files several wavelength bins have been co-added. 
This procedure has a twofold advantage because 
it allows us to (i) reconstruct the intrinsic pro- 
file and the corresponding SED and, (ii) to recog- 
nise whether a particular feature is resolved or not. 
The spectrum, extracted in a fixed width around 
the maximum intensity, is corrected by the ratio 
between the intrinsic and observed profile. Since 
for the LL2 segment the above procedure is gen- 
erally not as stable as for SL segments, we have 
preferred to fix Rq to the corresponding value de- 
rived in the nearby wavelength region of the SL 
segment. 

An estimate of the signal to noise (S/N) ra- 
tio was performed by considering two sources of 
noise: the instrumental plus background noise and 
the poissonian noise of the source. The former 



was evaluated by measuring the variance of pixel 
values in background-subtracted coadded images 
far from the source. The poissonian noise of the 
sources was estimated as the square root of the 
ratio between the variance of the number of e~ 
extracted per pixel in each exposure, and the num- 
ber of the exposures. The total noise was obtained 
by summing the two sources in quadrature and by 
multiplying by the square root of the extraction 
width in pixels. The corresponding S/N ratio at 
6/Lim is shown in column 6 of Table 1 . 

Finally we notice that the overall absolute pho- 
tometric uncertainty of IRS is 10%, while the slope 
deviation within a single segment (afl'ecting all 
spectra in the same way) is less than 3% (see the 
Spitzer Observer Manual). 

3. Results 

The final flux calibrated spectra of the selected 
Virgo cluster early type galaxies are shown in Fig- 
ures 1 and 2. 

3.1. Silicate emission from evolved stars 

In Figure 1 we have collected the thirteen galax- 
ies (76% of the sample) whose IRS spectra are 
characterised by the presence of a broad emis- 
sion feature around A r^lOfim that extends to- 
ward longer wavelengths. Thesc^ galaxies show 
neither PAH features nor emission lines. The 
observed spectra (solid lines) are superimposed 
on old SSP from Brcssan et al. (1998) normal- 
ized at A ~5.3/im. The dotted line is a 10 Gyr, 
Z=0.02 (solar metallicity) SSP computed with- 
out accounting for dusty circumstellar envelopes. 
Dashed lines from bottom to top refer to 10 Gyr 
SSPs with increasing metallicity Z=0.008, Z=0.02 
and Z=0.05, computed with dusty silicate circum- 
stellar envelopes. The models that account for 
dusty circumstellar envelopes show an extended 
feature due to silicate emission, which is very sim- 
ilar to that observed. The feature gets stronger 
at decreasing age and/or at incTcasing metallic- 
ity due to the corresponding higher dust mass loss 
rate of the SSP. Since, in addition to the match 
with the models, the analysis of the intrinsic spa- 
tial proflle indicates that the whole spectrum is ex- 
tended, we argue that the observed features are of 
stellar origin and most likely arise from dusty cir- 
cumstellar envelopes of mass-losing, evolved stars. 
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Fig. 1. — IRS spectra (solid lines) of passively 
evolving early-type galaxies in the Virgo cluster. 
Superimposed are SSP models from Bressan et al. 
(1998) normalized at 5.5/xm. The dotted line is a 
10 Gyr, Z=0.02 SSP computed without aceoimt- 
ing for dusty circumstellar envelopes. Dashed lines 
from bottom to top are 10 Gyr SSPs with metal- 
licity Z=0.008, Z=0.02 and Z=0.05 respectively, 
computed with dusty circumstellar envelopes. The 
dot-dashed line is a young (5 Gyr) metal poor 
(Z=0.008) SSP intended to show that also the 
MIR spectral region suffers from degeneracy. The 
bottom right panel compares the normalized con- 
tinuum subtracted spectrum of NGC 4365 with 
the normalized continuum subtracted spectrum of 
the mean outflow AGB star (dashed line, Molster 
et al. 2002) and the carbon rich star U Cam (dot- 
ted line, Sloan et al. 1998). For NGC 4365 the 
pseudo-continuum is a straight line that interpo- 
lates the observed spectrum between A=8//rn and 
A=13.5/xm. The spectrum of NGC 4339 is affected 
by poorly corrected droops caused by a star falling 
in the peakup. 



^ -2.30 \; ■, ^ 

-2.60 h '■-^ 

-2.70 j- , 1 

6 8 10 12 14 

X(/im) 

Fig. 3.— IRS spectra of NGC 4551 (dotted line) 

is compared with that of NGC 4365 (solid line). 
Spectra are normalized at A=5.3/um. The silicate 
features in NGC 4551 are less pronounced than in 
the case of NGC 4365, indicating a smaller dust 
mass-loss rate that can be attributed either to an 
older age or to a lower metallicity. 
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To corroborate this possibility wc compare, in Fig. 
1, the normahzed continuum-subtracted 10/Ltm sil- 
icate emission of the mean outBow oxygen-rich 
AGB star (Molstcr ct al. 2002) and of U Cam 
(a carbon rich star with SiC emission, Sloan et 
al. 1998) with that of NGC 4365. Though 
early type galaxies arc expected to harbour car- 
bon stars, given the wide metallicity spread within 
a galaxy, it seems that the dominant contribution 
comes from evolved M giants. More detailed mod- 
els, fully accounting for the expected mixture of 
evolved stars, will be presented in a forthcoming 
paper. 

The MIR view of early type galaxies proves to 
be a strong diagnostic for the population content 
of these galaxies. Recently Temi et al. (2005) no- 
ticed that the IRAC flux ratios at Sfxm and 3.6/xm 
or MIPS 24/im to IRAC 3.6/xm flux ratio remain 
fairly constant in early type galaxies which other- 
wise show different H/3 strength. They conclude 
that this disagreement supports a small rejuvena- 
tion episode. Although this is one of the possi- 
bilities invoked by Bressan et al. (1996) to ex- 
plain early type galaxies with strong Balmer line 
absorptions, caution must be paid before drawing 
definite conclusions. Indeed, we show in Figure 1 
that a young (5 Gyr) more metal poor (Z=0.008) 
SSP, dot-dashed line, is very similar to an old, 
more metal rich one. Thus even the mid infrared 
spectral region is degenerate and in order to break 
the age-metallicity degeneracy in passively evolved 
systems a careful combined optical (including pos- 
sibly NIR) and MIR analysis is required. To fur- 
ther illustrate the strength of this kind of analy- 
sis, we show in Figure 3 a comparison of the IRS 
spectra of NGC 4551 and NGC 4365. A recent 
optical spectroscopic study (Yamada ct al. 2006) 
indicates that NGC 4551 is significantly younger 
and more metal rich than NGC 4365. In this case 
we would expect NGC 4551 to be richer in bright 
mass-losing AGB stars than NGC 4365, and its sil- 
icate features to be more prominent. However, the 
opposite is observed, suggesting that NGC 4551 is 
either older or more metal poor (or both) than 
NGC 4365. Evidently, the effects of degeneracy in 
the optical can be strong (see e.g. Denicolo et al. 
2005, Annibah et al. 2006). 

We finally notice that NGC 4473 was observed 
by ISO (Xilouris et al. 2004) and shows spatially 
extended emission at 6.7 and 15/xm. These authors 



measured a 15/im excess with respect to SSP mod- 
els without dusty circumstellar envelopes. The ex- 
cess was interpreted as due to hot diffuse inter- 
stellar dust. IRS spectra, such as those presented 
here, permit the disentangling of the contribution 
of evolved AGB stars and the presence of inter- 
stellar dust. 

3.2. Active galaxies 

The remaining four galaxies (24% of the sam- 
ple) display different signatures of activity in the 
MIR spectra (Figure 2). These galaxies are clas- 
sified as active from optical studies (from AGN 
to transition Lincr-HII) at odds with the for- 
mer group. The spectra of NGC 4636 and NGC 
4486 (M 87) show emission lines ([ArII]7/um, 
[NeII]12.8Aim, [NcIII]15.5/im and [SIII]18.7/im) 
possibly of non-stellar origin. The broad contin- 
uum feature at lO/itm in NGC 4486 is not spatially 
extended and likely due to silicate emission from 
the dusty torus (Siebenmorgen et al. 2005; Hao 
et al. 2005). Line emission in NGC 4636 falls on 
top of the circumstellar emission SED and, as for 
NGC 4473, its excess at 15/xm is of stellar origin 
and not due to emission by hot diffuse dust as 
suggested by Ferrari et al. (2002). 

The spectrum of NGC 4550 shows PAH emis- 
sions features (at 6.2, 7.7, 8.6, 11.3 and 12.7 ^m) 
and the H2 S(5) 6.9/im and S(3) 9.66/im emis- 
sion lines. NGC 4435 shows a typical star-forming 
spectrum. A preliminary interpretation suggests 
that an unresolved starburst is dominating the 
MIR emission (Panuzzo et al. in preparation). 

4. Conclusions 

We presented Spitzer MIR IRS spectra of early 
type galaxies selected along the colour-magnitude 
relation of the Virgo cluster. 

We have reconstructed the intrinsic SED of 
these galaxies from the observed spatial profile 
sampled by the slits, via a careful analysis of PSF 
effects. In this way we are also able to differentiate 
between spatially resolved and unresolved regions 
within the spectrum. This provides independent 
support for the interpretation of their nature. 

Most of the galaxies (76%) show an excess at 
10/Ltm and longward which appears spatially ex- 
tended and is likely due to silicate emission. This 
class of spectra do not show any other emission 
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features. We argue that the 10/im excess arises 
from mass-losing evolved stars, as predicted by ad- 
equate SSP models. A detailed modelling of these 
features together with the analysis of combined 
optical, NIR and MIR spectra will be presented in 
a forthcoming paper. 

In the remaining smaller fraction (24%) we de- 
tect signatures of activity at different levels. We 
observe line emission superimposed on the stel- 
lar silicate features in NGC 4636, unresolved line 
and silicate emission in M 87 that likely originate 
in the dusty torus and unresolved PAH emission 
in NGC 4550 and NGC 4435. The latter galaxy 
displays the main characteristics of a nuclear star- 
burst (Panuzzo et al. in preparation). 

If we exclude M 87, which is a well known 
AGN, only two out of 16 early-type galaxies ob- 
served show PAHs, which corresponds to quite a 
low fraction (~12%) of the observed sample. It is 
premature to conclude that such a low fraction of 
galaxies with PAHs is representative of the cluster 
early-type galaxy population, especially if we con- 
sider that our investigation is limited to the bright- 
est cluster members (the upper two magnitudes of 
the colour-magnitude relation). A detailed com- 
parison of our results with those obtained for field 
galaxies will cast light on the role of environment 
in the galaxy evolution process. 

This work is based on observations made with 
the Spitzer Space Telescope, which is operated by 
the JPL, Caltcch imder a contract with NASA. We 
thank J.D.T. Smith for helpful suggestions on the 
IRS flux calibration procedure and the anonymous 
referee for useful suggestions. A. B., G.L. G. and 
L. S. thank INAOE for warm hospitality. 
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